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Abstract 

The nucleoli's tensor charges (isovector = 6u — 6d and isoscalar = 6u + 6d) are 
calculated using the QCD sum rules in the presence of an external tensor field. In addi- 
tion to the standard quark and gluon condensates, new condensates described by vacuum 
susceptibilities are induced by the external field. The latter contributions to g^ and g^ are 
estimated to be small. After deriving some simplifying formulas, a detailed sum rule analysis 
yields g^ = 1.29 ± 0.51 and = 1.37 ± 0.55, or 6u = 1.33 ± 0.53 and 6d = 0.04 ± 0.02 at 
the scale of 1 GeV^. 



PACS number(s):12.38Lg, 14.20Dh, 11.55Hx 



I. INTRODUCTION 



Studying nucleon charges such as the baryon charge, the axial charge, and the tensor 
charge from the underlying theory of the strong interactions, quantum chromodynamics 
(QCD), is obviously important to understand the properties of the nucleon as well as the 
QCD at low energies. On one hand, these nucleon charges, defined as matrix elements 
of bilinear fields {ip'j'^ilj for baryon charge, ip'y'^'y^ip for axial charge, and ipa'^'^ilj for tensor 
charge), are fundamental observables characterizing properties of the nucleon. On the other 
hand, these nucleon charges are connected through deep-inelastic sum rules to the leading- 
twist nucleon structure functions — -the spin averaged structure function fi{x), the quark 
helicity distribution and the transvcrsity distribution hi{x) [2, 3]. For example, 

the first moment of gi{x) is connected to the axial charge by the Bjorken sum rule [1]. 
Recently, Jaffe and Ji [3] showed that the first moment of hi{x) is related to the tensor 
charge Sip{ip — u, d...) of the nucleon through 



Thus, the tensor charge, like the axial charge, is one of the important nucleon observables. 

So far, there is yet no experimental data on hi{x) and on the tensor charge because hi{x) 
does not manifest itself in inclusive deep inelastic process and there are no fundamental 
probes which couple directly to the tensor current. A first measurement of h\{x) may be 
performed at HERA with semi-inclusive deep-inelastic scattering and at RHIC with trans- 
versely polarized Drell- Yan process [4] , from which one may extract the value of the tensor 
charge. Therefore, theoretical study of the tensor charge of the nucleon in QCD is very 
important. 

In Ref . [5] , we sought to understand the physical significance of the tensor charge and 
made estimates in the MIT bag model and the QCD sum rule approach. The QCD sum 
rule approach [6] has been used extensively in the past as a useful tool to study hadron 
masses, coupling constants, form factors, etc. (see [7] for a review). To evaluate the tensor 
charge and other charges, one makes an extension of the standard QCD sum rules [6] to 
avoid infrared divergences arising from zero-momentum transfer in bilinear quark operators. 
There are several equivalent formulations of the extension: the two-point correlation function 
in an external field, i.e., the external field QCD sum rule approach [8], and the three-point 
function approach [9]. In Ref. [5], we used the three-point function approach to calculate the 
tensor charge, including operators up to dimension 6. 

In the present paper we refine and extend our work in Ref. [5], in particular, we now 
include operators of dimension 8. To simply the calculation, we use instead the QCD sum 
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rule approach in external field [8]. The basic idea is this: Consider propagation of the 
nucleon current in the presence of an external tensor field Z^^. The two-point correlation 
function of the nucleon current is calculated up to terms linearly proportional to Z^y^ and 
the coefficient of which is the linear response of the correlation function to the external 
Zy^y field. The linear response is first calculated in terms of vacuum properties through an 
operator product expansion(OPE) [10]. In addition to the standard vacuum condensates [6], 
there are new condensates induced by the field Z^^^y. The latter represents the response 
of the QCD vacuum to the external field and can be described by vacuum susceptibilities. 
The susceptibility terms are exactly equivalent to the bi-local contributions in the three-point 
correlation function approach [5]. On the other hand, the linear response can be expressed in 
terms of properties of physical intermediate hadron states, where the tensor charge enters. 
By matching the two calculations in the certain momentum range, we extract the tensor 
charge. 

We first look at the sum rules in which there are no contributions of vacuum susceptibil- 
ities. By eliminating the nucleon coupling constant through the sum rule for mass|TT[], we 
obtain the simplifying formulas of the tensor charges. 



= ^ (2) 
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at the scale /x^ = M^. [The leading-log evolution of the tensor charge is 



^?^(/i^) = [«(^^)/a(/i^)]"/^^(;.^) 

where indices a = v,s,u,d specify = 6u, = 6d, g^ = 6u — 6d, and g^ = 6u + 6d, 
and Uf is the flavor number.] Substituting {qq) = -(240MeV)^, {glC^) = 0.47GeV'^,mg = 
-{qgccr ■ Gq)/{qq) ^ O.SGeV^ and mjv = 0.94GeV, we get 5u = 1.29, Sd ^ 0.02, g^ = 1.27 
and g^ = 1.31. The uncertainty from the vacuum condensates alone is at the level of 20%. 
In order to extract these tensor charges with a better accuracy, we analyze the sum rules in 
a more careful way. The result is consistent with the simplifying analysis presented above. 

The paper is organized as follows. In Sec. II we derive the nucleon sum rules in the pres- 
ence of an external tensor field. The sum rules are analysed in Sec. Ill, where we first derive 
the simplifying formulas for the tensor charges, then estimate the vacuum susceptibilities, 
and finally make a detailed sum rule analysis. The summary and conclusions are given in 
the last section. 
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II. QCD Sum Rules for the Nucleon in An External Tensor Field 



To derive the QCD sum rules, we consider the two-point correlation function in the 
presence of an external constant tensor field Z^y, 

U{Z,,,p) = i J d'xe'P^{0\T{ri{x)fi{Omz: (3) 
where rj is the interpolating current for a proton [11], 

= eabc[ul{x)C-fu,Ub{x)]-f5-f^'dc{x) , (4) 

{0\vmNip)) = \NVNip) , (5) 

where a, b and c are color indices, superscript T means transpose, C is the usual charge- 
conjugation matrix, and vn{p) the nucleon spinor normalized to vv = 2mjv. To construct 

the sum rule for tensor charges, we expand the correlation function (3) to the first order in 
7 

U{Z^,,p) = Uo{p) + Z^^U^^'ip) + (6) 

where IIq is the correlation function which can be used to construct the nucleon-mass sum 
rules. In the following, we focus on the linear response function n^'^(p). 

The QCD sum rule is based on the principle of duality, which postulates the correspon- 
dence between a description of correlation functions in terms of hadronic degrees of freedom 
(phenomenological side) and that in quark and gluon degrees of freedom (QCD theory side). 
The phenomenological representation of the correlation function in an external Z^i^ field can 
be expressed in terms of physical intermediate states as 

U(Z,,,p) = g^Z.^lpa'^'^p + m\a^^ + m^{p, a^^}\ + • • • , (7) 

(p^ - m^) 

where the nucleon 's tensor charge is defined by 

(7V(p, s)| j;,(0)|7V(p, s)) = g^v^{p, s)a^,v^{p, s) , (8) 
where J^^ = 'ipaf^nv'^a with a = v,s,u, or d, specifying J'"l^^ = ua^yU =F dc^nud and J^^'^) — 



ua. 



fJ.U 



u{da^yd). 



The theoretical side of the correlation function (3) can be computed in the deep Euclidean 
region using OPE. Matching the two results obtained above, we can extract the tensor charges 
in terms of QCD parameters and vacuum condensates. 



A. Operator Product Expansion 
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We now compute the OPE for the nucleoli- current correlation function in an external 
tensor field Z^y. The coupling between quarks and the external antisymmetric tensor field 
Z^y is described by an additional term 

A£ = g^qa^'^qZ^y (9) 

in QCD Lagrangian, where the coupling constant depends on the quark type as well as 
the field Z^y{gq = gu — —ga for isovector type coupling, gq = gu — ga for isoscalar type). 
The equation of motion for quarks is now 

(zZ) + y,Z^,ang = , (10) 

where — dfj, + igcA^^^,gc is the QCD gauge coupling, A" are the Gell-Mann matrices, 
and current masses of up and down quarks are neglected. 

To calculate the Wilson coefficients in OPE [10], we construct the coordinate-space quark 
propagator in the presence of the external field Z^j^^. Following the method of Refs.[8] and 
[12], we find 

^? = (0|rg«(x)g-5(0)|0) 

10 yx)^J alx'^Z (a'"') --4:X^xZ (a"")-] 



+ ^9c—G';Jxa'"' + a^^x,,, 

1 5' 



--5^^5,Mm - -^Z,ygqx{qq){^ni3 

Xab X 

-i—g^l^qq)Z^y{xa^^ + a^^x),^ + ^5,ix\Q\qg,a ■ Gq\Q) 

Sjab 

+—g^{qq)Z^y[x\K - 20{an^J + 2a;%(«: + C){cr^%] + ■■■. (11) 

The first three terms in Eq.(ll) are the perturbative propagator for a massless free quark 
and its interactions with the external Zf^^ field and the vacuum gluon field G^^. The next 
five terms are non-perturbative, arising from the quark and the mixed quark-gluon conden- 
sates, and the condensates induced by the external field Z^,,. In fact, due to breakdown of 
Lorentz invariance, the vacuum expectation values (qcr^yq), {qgc^G'^^q) , (qgclbG^yq), where 
G^u = \^^jLursG'^^^^, become non-zero. To characterize these vacuum response, we define 
the induced susceptibilities %, and C, through 



(0|ga^^g|0)z = gqXZ^ii.{qq) , 
{^\Q9c^GIu(1\^) z = gqKZ^y{qq) , 

{^\qgclbG^yq\^)z = -igqCZi,y{qq) . (12) 
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Using 

(0|r(r^(x)7^(0))|0) = -26»''^6«'''''='rr{5(x)t'''7.^^'5(^)r'C7M}757'^'5(^)f 7'^75, (13) 

and Eq. (11), it is straightforward to carry out the OPE. The corresponding Feynman 
diagrams are shown in Figs. 1 and 2. The results for the correlation function (3), including 
the contributions up to dimension 8, is 

n(Z^,, p) = Z^, [Wipa^^p + W2a^'' + W^{p, a^'^} + • • • ] (14) 



where 



^'^{m? -^{(l<l){(19cO-Gq) , (15) 



^ J^pW) + ^(4«: + C)(?-?)M-/) 
Ws = |^(gg)/n(-/) + 24^^^(-//A')(g^ca-Gg) 

^''-,x{mf + ■wrix{m){(igc(T-Gq) 



B. Sum Rules 

The QCD sum rules are obtained by equating the OPE results in Eqs.(14)-(17) and the 
phenomenological description in Eq.(7), and performing Borel transformation. The renor- 
malization scale dependence is also included properly. There are three different Dirac struc- 
tures: chiral-odd Z^j^i^pa^^p and Z^^a^^, and chiral-even Z^i,{p,a^^}, each of which can be 
used to construct a sum rule and extract the tensor charges. Note that the phenomenological 
representation of Il{Zf^i,,p) in Eq.(7) contains the double-pole as well as various single-pole 
contributions. The latter terms can be treated in several ways. Here we eliminate them by 
multiplying both sides of the sum rules by — p^ before Borel transformation. We find 
the sum rule for the structure Zn^pa^^p, 

M\mlE, - M^E,) + ^m^LV- - ^(1 + ^) 



24LV9 y ^ ^ 16 ^ 576LV9^ 

I 9uO' /-. ^AT w 16/27 _ Sgf^mpg^ m% 4/27 

2 ^ M^' 8M2 2M^' 

= ^?(/^o)/?l^e-"^-/^^ (18) 
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for 



^bM'{m%Eo-M'E,)+'^a'ml 
9o D 



(19) 



and for Z^^{p, a^"}, 



'-^M\m%E, - M-E,)L^'-' + '-^^^{Hml/ K^) - 1) 





(20) 



Here a = -{27rf{qq), b = {g^G^}, mla = {27rf{qg,a ■ Gq), l3% = (27r)^A^/4, and M is 



the contributions from excited states, where sq is an effective continuum threshold. The 
anomalous dimension of various operators, including that of the tensor current, is taken 
into accout through the factor L = ln{M'^ / h?) /ln{^^/ h?') [6, 11], where /xq is the inital 
renormalization scale (500 MeV) and A = Aqcd is the QCD scale parameter (100 MeV). 
The tensor charges in Eqs.(18)-(20) are renormalized at the scale /iq. 



Let us now analyse the sum rules derived in the previous sections and extract the tensor 
charges g^, g^ and Su, 5d. In principle, one expects to obtain the same result from each of 
the sum rules in Eqs.(18)-(20). In practice, however, some sum rules work better than the 
others, as was evident in similar studies [8, 14, 15]. In our case, the sum rule in Eq. (18) 
from the structure pa^^p is preferred for the following reasons. Firstly, from Eqs.(14)-(17) 
the structure pa^^p contains extra powers of momentum in the numerator compared with 
structures a^^ and {p, a^^}. Hence, the resulting sum rule converges better after Borel trans- 
formation. Secondly, the sum rule in Eq. (18) is known better phenomenologically, because 
the terms involving the poorly-known susceptibilities are numerically small compared with 
the dominant dimension-6 and 8 terms. Hereafter, we discard the sum rule in Eq. (19). 

A. Simple Formulas for the Tensor Charges 



the Borel mass. The factors K 



n 



1 ~ S~fe ^ with X = sq/M"^ account for the sum of 



III. Sum Rule Analysis 
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Consider the sum rules (18) for the structure pa^^p. As a first approximation, we dis- 
regard the anomalous dimensions, let M = rriN, and use /S^e^^Af/^-'^^ |^^2=j„2 = — 7r^m^(-u-u) 
from the simplified mass sum rules [11]. In the case of gu = l,gd = 0, we find a simple 
formula for the tensor charge Su, 



4(27r)2(gg) 9ml ^ 



Similarly, from the sum rule (20) with the structure {p, cr'^''}, we get 

^d-^-^. (22) 

The result obtained above is independent of the susceptibilities x, k and Thus, one 
can get a simple estimate of the tensor charges without knowing them. Substituting the 
vacuum condensates into the above equations, we find Su 1.29, Sd 0.02, or 1.27 
and gj, ~ 1.31 (at /x^ = "^at)- Interestingly, 5d is much smaller than 5u, or isovector gj, and 
isoscalar g^. tensor charges have similar sizes. 

B. Estimates of Susceptibilities 

To make a more detailed analysis, we need the values of the susceptibilities x, k, and 
C defined in Eq.(12). These susceptibilities correspond to the bi-local contributions in the 
three-point correlation function approach [5, 9]. In fact, using Eq. (9), it is easy to show 

{0\qa^''q\0)z = g,Z^, d''xe'^%Q\T{qa^pq{x), qa^^^qM . (23) 

Combining with Eq.(12), we obtain 

X(gg)= lim^/ciSe'f-(0|r(ga,,^g(x),g(7"^g)|0)=^n^(0) , (24) 

where II^(O) is a bi-local correlator defined in the three-point correlation function approach 
(see Eq.(23) of Ref.[5]). Similarly, we have 

<qq) = \j d'x{0\T{qg,^GZ^q, ga«^)|0)=^n.(0) , (25) 

cm = -lJd'x{0\T{qg,^,G^pq, g(7°^)|0) = -^0^(0) . (26) 

Substituting the above expressions for x, and ( into Eqs.(15)-(17) and (18)-(20), we obtain 
the same result obtained in the three-point correlation function approach [5]. 

Eq.(24)-(26) provide useful formulas to determine the susceptibilities Xi and (. To find 
X, we assume p(l ) and 5(1+^) meson dominance in Eq. (23). Estimating the coupling 
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constants of relevant currents with the meson states again using the QCD sum rules, we get 
[5] n^(0) ^ (0.15GeV)2, or 

Xa = -^-^n^(0) ^ -O.lSGeV^ (27) 

As pointed out in Ref.[5], the small 11,^(0) comes from the cancellation of the two resonances. 
However, even if p meson does not couple to the current qa'^^q [13], we have 11-^(0) ~ 
(0.21GeV)2 and xa ~ -0.29GeV^ which are still small. 

The susceptibilities k and C, can be estimated in a similar way. To get just an order of 
magnitude, we assume only B-meson couple to the current qa^^q. We find 

n,(0) = -n.(0) ^ e<l^\^-^{mqq) + ^M^')) , (28) 

where M and are the Borel mass and the vector meson mass, respectively. Substituting 
{mqq) ^ - {100MeV)\ (g^^G^) ^ 0.47GeV^m^ ^ 1.2GeV,M ^ IGeV, we get n^(0) = 
-n«(0) « 0.017GeV^ or 

Ca = KOfti O.lOGeV^ . (29) 

The above estimates of k, and ( are admittedly crude. Fortunately, as we shall see from 
the anabelow, the effects of susceptibilities on the tensor charges and g^ are small. 



C. Sum Rule Analysis 

Now we perform a standard sum rule analysis for Eqs. (18) and (20). Multiplying 
both sides of equations by we study g^{p) as function of Borel mass M^. In our 

calculation, sq is taken to be 2.3GeV^, (3% = 0.26GeV®, and other parameters are the same 
as those following Eq. (2). The solutions for 5u{iJ?) and 5d{jj,'^) are plotted as function of 
in Fig. 3. Su from the two sum rules are shown in solid and dashed hues, Sd are shown 
in dot-dashed and short-dashed lines. We take the spread of the two sum rule predictions 
as our theoretical uncertainty. Averaging the two results, we have 

5m = 1.33 ±0.53 , (5ci = 0.04 ± 0.02 (30) 

at the scale /j,^ — — (IGeV)^. Correspondingly, 

= 1.29 ±0.51 , = 1.37 ±0.55, (31) 

which are consistent with the simplified results early. 
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IV. Summary and Conclusions 



In this paper we studied the nucleon's tensor charges (isovector one g^, isoscalar one g^, 
and 6u and 6d) by means of QCD sum rule approach. The calculations include terms up to 
dimension 8, which presumably gives more reliable predictions than our previous result. Our 
final numbers for the proton's tensor charges are g^ = 1.29 ± 0.51, g'^ = 1.37 ± 0.51, 6u = 
1.33 ± 0.53, 6d = 0.04 ± 0.02 at the scale of 1 GeV^. 

In obtaining the above results, we used different methods to analyze the sum rules for 
g^. The first one is to derive the simplifying formulas for the tensor charges at /i^ = m^. 
The second method is to analyze the sum rules with different chiral structure, from which 
we obtained the tensor charge at the scale of 1 GeV^. Both analyses seem to be consistent. 

Our QCD sum rule calculation shows that 6u ^ 6d, which means that the up quark 
dominates the contribution in a transversely-polarized proton. Furthemore, the isovector 
tensor charge g^ of the proton is comparable in magnitude with the isovector axial charge 
gA- To the contrary, the isoscalar tensor charge is markedly different from the isoscalar axial 
charge. We hope that the future experimental measurement will test the present prediction. 

At last, we wish to point out that since Sd is small, the precise value of which, including 
its sign, is difficult to calculate. Futhermore, the sum rules presented here have no stable 
solutions, which may refiect that some relevant physics is missing in the OPE. For instance, 
instanton physics, important in the case of the nucleon mass ||16|, might stablize these sum 
rules and improve the present prediction. 

We thank Xuemin Jin for useful discussions. H.X. He was supported in part by the 
National Natural Science Foundation of China and the Nuclear Science Foundation of China. 
X. Ji was supported in part by funds provided by the U.S. Department of Energy (DOE) 
under cooperative agreement No.DF-FC02-94ER40818. 
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Figure Captions 

Fig.l. Diagrams for the calculation of the Wilson coefficients of the correlation function 
correspnding to the structures pa^^p and a^'' . 

Fig. 2. Diagrams for the calculation of the Wilson coefficients of the correlation function 
corresponding to the structure {p^a^"}. 

Fig. 3. The solutions of tensor charges from Eqs. (18) and (20). The solid and dashed 
lines represent 6u, and the dot-dashed and short-dashed lines represent 6d, with upper 
curves from Eq. (18). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9607408vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9607408vl 



